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ABSTRACT 

The  finite  difference  (FD)  technique  is  used  to  investigate  the  air  gap  effect  on  the  coupling 
between  two  planar  microstiip  lines.  The  strip  lines  are  embedded  between  a  dielectric  substrate 
and  a  dielectric  overlay.  The  air  gap  between  the  substrate  and  the  overlay  is  present  due  to  the 
thickness  of  the  strips  and/or  the  roughness  of  the  planar  surfaces.  A  quasi-static  FD  technique 
based  on  the  solution  of  Laplace's  equation  is  used  to  solve  for  the  potential  distribution  in  the 
region  around  the  strips.  The  FD  mesh  is  truncated  by  an  artificial  boundary  to  allow  for  a 
numerical  solution  with  the  available  computer  resources.  At  this  artificial  boundary  an 
approximate  boundary  condition  is  used  to  truncate  the  FD  mesh.  From  the  potential  distribution, 
the  charge  and  the  capacitance  and  inductance  matrices  of  the  stripline  system  are  evaluated.  The 
coupling  or  cross-talk  between  the  lines  are  then  defined  and  computed.  The  effects  of  the 
conducting  strips  thickness  and  the  properties  of  the  materials  between  the  strips  on  the  coupling 
are  illustrated  through  numerical  examples. 
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INTRODUCTION 
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Over  the  past  few  decades,  microstrip  transmission  lines  have  been  analyzed 
extensively.  Reducing  the  coupling  has  become  very  important,  since  as  the  clock  rates  in 
computers  increase  so  does  the  effect  of  coupling  or  cross  talk  between  the  lines.  The  main 
purpose  of  this  work  is  to  analyze  the  effect  of  the  air  gap  on  the  coupling  between  two 
planar  microstrip  lines.  Two  different  geometries  are  analyzed  using  the  finite  difference 
(FD)  technique.  Figures  1  and  2  show  the  geometries  where  the  striplines  are  perfectly 
embedded  between  the  substrate  and  the  overlay  in  contrast  with  Figs.  3  and  4  where  an 
air  gap  is  present  due  to  the  thickness  of  the  strips  and/or  the  roughness  of  the  planar 
surfaces.  Following  the  approach  presented  in  [1],  a  FD  technique  based  on  the  solution 
of  Laplace's  equation  is  used  to  find  the  potential  distribution  around  the  strips.  From  the 
potential  distribution  the  charges  on  the  strips  are  computed  using  Gauss'  law  and  then  the 
capacitance  and  inductances  matrices  are  evaluated.  The  capacitive  coupling  between  the 
lines  is  then  computed.  Numerical  data  are  obtained  for  these  strip  line  geometries  with 
and  without  air  gap  in  order  to  show  the  effects  of  the  air  gap  on  the  coupling  between  the 
lines.  The  developed  program  is  capable  of  solving  different  geometries  of  multi 
conductor  transmission  lines  arbitrarily  positioned  in  an  inhomogeneous  dielectric  space. 
A  FD  mesh  generator  program  is  also  developed  to  act  as  a  preprocessor  for  the  FD 
routine.  Lagrangian  polynomials  are  used  to  approximate  the  second  order 
differential  equation  in  Laplace's  equation. Then,  the  integral  form  of  Gauss's  law  is  used 
to  derive  a  general  expression  for  the  scalar  potential  at  every  point  on  the  FD  mesh.  An 
Artificial  Boundary  Condition  (ABC)  is  used  to  truncate  the  finite  difference  mesh  to  allow 
for  a  numerical  solution  with  the  available  computer  resources.  The  resulting  expressions 
for  the  potential  at  the  FD  mesh  points  are  cast  into  a  matrix  equation.  Special  arrangement 
of  the  points  allow  the  use  of  a  banded  matrix  solver  for  the  solution  of  the  related  system 
of  equations.  From  this  solution  of  the  matrix  equation  (i.e.  the  potential  distribution),  the 
electric  field  vectors  are  computed. 
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it  rate 


1.  A  transmision  line  with  two  strips  between  the 
substrate  and  extended  overlay  (no  air  gap). 


.  2.  A  transmision  line  with  two  strips  between  the 
substrate  and  the  dielectric  overlay  (no  air  gap). 


Artiticial  Boundary 


Fig.  3.  A  transmision  line  with  two  strips  between  the 
substrate  and  extended  overlay  (with  air  gap). 


_ Artificial  Boundary 


Fig.  4.  A  transmision  line  with  two  strips  between  the 

substrate  and  the  dielectric  overlay  (with  air  gap). 


CHAPTER  n 

FINITE  DIFFERENCE  FORMULATION 
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A)  Finite  Difference  Approximation 

Through  the  years,  many  methods  have  been  developed  for  the  analysis  of 
microstrip  transmission  lines.  Among  those  is  the  widely  known  finite  difference 
technique.  For  the  quasi-static  2D  problems  as  shown  in  Figs.  1  to  4,  the  potential  V  can 
be  described  using  Laplace's  equation  by: 


-  +  -  =  u 

dx^  dy^ 

The  finite  difference  technique  approximates  the  derivatives  involved  in  Laplace's 
equation.  This  numerical  approximation  is  based  upon  calculating  the  potential  at  a  certain 
node  as  a  function  of  neighboring  grid  points,  or  nodes. 

In  general,  assuming  three  consecutive  nodes  L,C,R,  as  shown  in  Fig.  5, 
Lagrange's  form  of  interpolating  polynomial  is  described  as  [2] 


f(x)‘‘PJx)- 


(X  X^)  (x  Xjj) 


(.X-Xj)  (X-Xj) 


(x-Xj)  (x-x^) 

(Xjj-X^)  (Xj^-X^) 


Ax^  (2) 


where  P2(x)  represent  a  second  degree  polynomial  which  coincides  with  the  exact  values 
of  the  function /fxj  at  the  three  nodes  L,C,R.  In  order  to  apply  Laplace's  equation  over 
the  potential  V(x,y),  the  second  order  derivatives  need  to  be  approximated  in  both  the  x  and 
y  directions.  The  Lagrange  interpolating  polynomial  will  then  have  to  be  differentiated 
twice,  that  is 


<y(x)  (X-X^,)  t  (X-Xj^)  ^  (X-X^)  t  (x-Xj^) 

dx  (x^-Xp)  (x^-Xjj)  (x^-Xj)  (x^-x^) 


(x  -  x^)  +  (x  -  x^) 
(x^-x^)  (x^-x^) 


Axj,) 


(3) 
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dx^  (X^-X^)  (X^-X^)  (Xp-X^)  (X^-Xjj)  (Xjj-X^  (Xjj-X^) 


Similar  expressions  for  the  derivatives  with  respect  to  y  can  be  obtained. 

The  boundary  conditions  for  the  microstrip  transmission  lines  shown  in  Figs.  1  to 
4  are  given  by 

V  =  Vj  on  strip  1 

V  =  V2  on  strip  2 

V  =  0  on  the  ground  plane 

£D-ds  =  0  on  the  dielectric  interfaces 


where  Vj  and  Vj  are  arbitrary  constant  potentials  and  s  is  the  closed  surface  around  a  grid 
point  on  the  interface. 

The  last  condition  is  based  on  Gauss's  law  for  the  electric  field,  i.e. 


jD-ds  -  jeE-ds-  -  0 


(5) 


where  is  set  to  zero  since  no  free  charges  exist  on  the  dielectric  boundaries.  The 
surface  integration  in  equation  (5)  can  be  replaced  by  a  contour  integration,  because  the 
microstrip  transmission  lines  considered  are  two-dimensional  structures,  and  the  solution 
of  the  potential  is  independent  of  z  [3]. 
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Substituting  E  =  -W  in  equation  (5)  yields  ; 


dV 


(6) 


where  dV/dn  denotes  the  derivative  of  V normal  to  the  contour  as  shown  in  Fig.  6. 
It  should  be  noted  here  that  the  normal  is  defined  as  the  vector  pointing  out  of  the  area 
enclosed  by  the  contour  I  Figure  6  shows  a  general  point  (C)  which  is  surrounded  by  four 
different  dielectric  materials,  this  is  the  most  general  case,  and  the  mesh  generator 
program  provides  the  data  about  all  these  dielectric  materials  around  such  point.  After 
applying  equation  (6)  to  the  contour  around  that  point,  an  expression  is  obtained  for  the 
computation  of  the  voltage  at  these  type  of  general  points,  that  is 


(^1*^2)  1 

+ 

yT~y  c ) 

Vr 

(^f^2)j 

2 

1 

1 

2 

(7) 


The  above  equation  deals  with  any  kind  of  point  because  the  permitivities  of  the  regions 
around  point  C  do  not  need  to  be  all  different.  If  all  dielectric  regions  are  the  same  then 
the  point  is  in  a  homogeneous  medium,  and  in  inhomogeneous  cases,  the  point  will  be 
located  at  a  surface  or  line  interface  between  two  or  three  different  media,  respectively. 
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B)  ArtiHcial  Boundary  Condition  at  the  outer  boundary 
For  the  transmission  lines  shown  in  Figs.  1  to  4,  approximate  boundary  conditions 
should  be  applied  at  the  outer  artificial  boundary  in  order  to  truncate  the  mesh  so  that 
numerical  analysis  could  be  carried  out.  The  equation  used  on  the  general  type  point  (type 
C)  can  not  be  used  here.  The  mesh  generator  program  defines  a  new  type  for  the  points 
on  the  boundary,  type  "B",  which  is  subdivided  in  subtypes  depending  on  the  location  of 
the  boundary  point.  At  these  points  the  approximate  boundary  condition  proposed  by 
Khebir  et.  al  is  used  [4].  The  series  expansion  for  the  electric  potential  is  defined  as  [5] 


F(p,(t))  =  C,.^,((i>)liip.5: 

n  1 


p” 


(8) 


where  the  constant  term,  Q,  is  dropped  because  the  potential  at  infinity  is  zero.  The 
second  term  in  equation  (8)  is  also  dropped  because  the  total  charge  enclosed  by  the  mesh 
at  the  artificial  boundary  is  also  zero.  Having  set  Q  and  Aq  to  zero,  equation  (8)  becomes 


^2(4>)  A0) 

F(p,<j))  =  - + - —  ♦ - —  + 

P  P'  p' 


(9) 


Now,  equation  (9)  is  differentiated  with  respect  to  p  ,  to  obtain: 
5F(p,4>)  3^3(4>) 

dp  p^  p^  p'* 


(10) 


Then,  multiply  equation  (9)  by  i/p  and  add  it  to  equation  (10): 


dV  I  y.  -^1  ■'^2  "^3  ^*^3 

’  7^"^ V ''  p"  p'  p' 


(11) 


Finally,  rearranging  equation  (1 1)  yields, 

^ 

5p  p  p*  p^  p*  p* 


(12) 
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As  an  approximation  to  equation  (12),  the  right  hand  side  is  assumed  to  be  zero.  This 
truncation  amounts  to  an  error  in  the  analysis  which  will  be  reduced  as  the  value  of  p  at 
the  outer  boundary  is  increased,  i.e.,  as  the  outer  boundary  of  the  mesh  is  moved  away 
from  the  strips. 

The  chain  rule  can  be  used  to  write  dVIdp  as: 

ap  '  dx  dp*  By  dp  ^  ^ 


where  p(x,y)  =  .  By  dropping  the  right  hand  side  of  equation  (12)  and 

substituting  the  result  into  equation  (13),  one  obtains 
V  dV  dx  dV  dy 

■  ,  3  *  (14) 

p  ax  ap  ay  ap 


Since  x  =  p  cos  (cj))  and  y  =  p  sin  ((j)),  the  derivatives  of  p  with  respect  to  x  and  y  are 


thus  defined  as: 


—  =  cos  ((|))  =  —  , 

ap  p 


Substituting  the  results  from  equation  (15)  into  equation  (14)  yields 


dV  dV 
-V  =  X  —  *y  — 

dx  dy 


Equation  (16)  can  be  rewritten  as 
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dv  ar'l 

—  =  —  V*y  — 
dx  X  \  dy  ^ 


(17) 


Equation  (17)  can  be  implemented  along  the  left  and  right  boundaries  of  the  mesh.  Figure 
7a  shows  the  grid  point  system  and  the  notations  for  boundary  points.  As  an  example,  on 
the  right  boundary  (type  "35"  points)  of  the  grid  point  Q  resides  on  the  boundary  and 
is  a  point  outside  the  outer  boundary.  Because  Rj  lies  outside  the  mesh,  it  is  not  possible 
to  enforce  Laplace's  equation  (1)  at  node  Q  in  the  usual  manner.  Instead,  equation  (17) 
is  used  to  determine  the  potential  at  R^,  i.e.  Vg,  in  terms  of  V^,  Vc,  Vj.,  and  Vg.  Then, 
once  Vg  is  known  equation  (1)  can  be  enforced.  The  normal  derivative  of  the  potential  at 
point  C5  on  the  right  boundary  is  defined  as 


dVc  Vg-V^ 
dx 


(18) 


where  2  A  =  x^^-x^  .  By  substituting  the  expression  derived  in  equation  (17)  into  equation 


(18),  one  obtains: 


r  =  F- 

R  L 


2L 


dy 


(19) 


The  term  dV^/dy  is  then  replaced  by  its  finite  difference  approximation  to  give 


2A 

r/  1/ 

V  ^ 

(yc->'r)  y  <yc-yB>  y  (yc-ys)  *  (yc-yi)  y  ^ 

^c*yc 

^  (>'s->'c)  ^  (yr->'B)  Or->'c)  ^  (yc-yB>  (>'c->'r) 

(20) 


Left  Boundary  (B1  type) 
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/ 

/ 
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Top  Boundary  (B3  type) 
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Right  Corner  (B4) 


*T, 


-•-Lv-®e,--*-Rv . 


B, 
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rx. 


Lj  fCj  •Rs 

Yv 


•L, 
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O) 

a 

u 

CO 

-O 
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O 

« 
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’  Bottom  Bottom  Boundary  (B7  type) 
Left  Corner  (B8) 


tT, 

^tr®C5  •Re 

•B 

Bfi  I 

Bottom 

Right  Corner  (B6) 


Fig.  7a.  Artificial  boundaries  showing  the  points  that  lie  outside 
of  the  FD  mesh  for  each  type  of  boundary  and  comer. 
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which  can  be  rearranged  to  form 


,  =  F,  +  -  1  +  w - 

[  *c  [  (yc-yB>  C>'c->'r)j. 

,V  ^c-y^  \ 

V  ^B-yc>(yB-yT^)\  yc  [  (yT-y^(yT-yc)j _ 

For  simplicity,  equation  (21)  can  be  written  as: 


(21) 


c,v,.  (22) 

where  Q  ,  ,  Q  ,  and  Cj-  represent  the  coefficients  of  the  potentials  as  defined  in 

equation  (21).  Now  that  Vg  is  expressed  in  terms  of  V^,  Vc,  Vg  ,  and  Vj- ,  Laplace's 
equation,  (1),  can  be  used  to  relate  the  potential  at  the  grid  points  on  the  boundary  to  their 
surrounding  nodes.  Applying  equation  (1)  at  node  Q  gives 

■  0  (23) 

where 
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Substituting  the  expression  for  in  (22)  in  to  (23)  yields 

.  V^{K^C^*K^-\*V,[K^C,.K,]  -  0  (25) 

This  is  the  equation  that  is  enforced  along  the  right  side  of  the  outer  boundary.  The  method 
used  along  the  left,  top  and  bottom  sides  of  the  outer  boundary  is  similar,  as  shown  in  Fig. 
7a.  Thus,  on  the  left  boundary  (type  "Bl"  points)  we  have 

r^[K^C,.K^].V,lK^C^,K^-\.V,[K^C,.K,].r^[K^C,.K,]  .  0  (26) 

where 

2A  { ,  (yc-ys)  *  (yc-yr) 

-  1  +  - 

V  (yc-ys)  ^c-yj) 

Oc->^r)  ] 

C  ] 

^  [  >'c  i  Cyr->'c)  j 


and  Kl,  Kji,  Kc,  Kj,  and  are  defined  in  equation  (24). 

In  order  to  apply  the  same  numerical  procedure  on  the  top  side  of  the  outer  bdundary , 
equation  (17)  should  be  rewritten  as 


dv  if-,  ar"! 

dy  y[  dx) 


(28) 


Starting  from  equation  (28),  the  same  analysis,  as  for  the  left  side,  is  repeated  to  obtain 
an  expression  for  the  top  side  (type  "B3"  points)  as 


and  again  Kj^,  Kc,  Kj,  and  are  defined  in  equation  (24).  From  equation  (28)  and 
following  the  same  procedure,  an  expression  is  obtained  for  the  lower  boundary  ("B7"  type 
points) 


F  [KC,  .  KA  .  [K^C^  .  K^] .  [K^C^  .  KJ  .  F,  [K^C,  ^K^]-0 


where 
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1  +x^ - 

^  (^C-^L^  ) 

^  (x^-x^)(x^-x^)j 

^  (^C-Xjj)  "l 

^  (x^-x^)(x^-x^)j 


(32) 


where  A  is  defined  as  A=(jj.-j'c)- Again  these  coefficients  must  be  substituted  in  the 
general  equation  to  obtain  an  expression  for  the  points  on  the  lower  boundary.  At  the  four 
comers  of  the  artificial  boundary  (types  B2,  B4,  B6,  and  B8),  another  procedure  is  used. 
Figure  7b  shows  a  sketch  of  these  comers  and  the  neighboring  nodes.  For  example,  at 
right  upper  comer,  the  potential  at  C  will  be  related  to  the  potential  at  the  left,  and  at 
the  bottom,  V^.  Because  equation  (12)  is  based  on  the  first  order  expansion  of  equation  (9), 
the  potential  at  the  comers  will  also  be  related  to  only  7/p.  Thus  =  p  Vfp,4>J  ■ 

Using  linear  interpolation,  the  value  of  Aj  at  can  be  related  to  its  value  at  4>b  and  ^  as 


(34) 
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C4 


B4 


C6 


Fig.  7b.  Artificial  boundaries  showing  points  at  the  comers 
and  points  on  the  boundary  from  which  the  potential 
at  Cn  may  be  computed  by  interpolation. 
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-Pl 


(39) 


C)  Construction  of  the  Matrix  Equation. 

The  equations  derived  above,  describe  the  potential  at  a  node  by  relating  it  to  the  potential 
at  neighbor  nodes  in  the  mesh.  Therefore,  application  of  those  equations  at  all  the  nodes 
in  the  mesh  leads  to  a  set  of  simultaneous  equations  that  can  be  described  as  the  matrix 
equation  [ri][X]  =  [fi].  This  matrix  equation  is  like  the  one  found  in  [6].  The  Matrix  [ri] 
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is  diagonally  dominant  with  the  elements  away  from  the  diagonal  being  zeros.  Therefore, 
[4]  is  a  banded  matrix,  and,  like  in  [6],  it  is  desirable  to  take  advantage  of  this  fact  when 
solving  the  system  of  equation.  Linpack  Subroutines  [7],  SGBCO  and  SGBSL  are  used  to 
solve  for  the  potential  distribution  on  the  mesh,  the  column  vector  [X].  Once  the  potential 
is  known  for  every  node  on  the  FD  mesh  then  the  characteristic  parameters  of  the 
transmission  line  can  be  determined. 


CHAPTER  m 

COMPUTATION  OF  THE  TRANSMISSION  LINE  PARAMETERS 
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A)  Introduction 

After  solving  Laplace's  equation,  the  potential  distribution  on  the  FD  mesh  is 
known.  From  the  potential  distribution  the  electric  fields  on  the  mesh  may  be  computed. 
The  characteristics  of  the  transmission  line  may  also  be  computed  from  the  potential 
distribution  [6]. 


B)  Computing  the  Electric  Field  Vectors  at  the  Mesh  Points 
Once  the  potential  distribution  on  the  mesh  has  been  computed,  the  field 
distribution  inside  the  FD  mesh  domain  can  be  found.  The  derivation  starts  from  the 
following  definition 

E=-W  (40) 


where  at  any  node  we  get 


E  =  - 


X 


(41) 


then  the  finite  difference  approximation  given  by  (2)  find  (3)  is  used  to  obtain  the  electric 
field  components  at  any  mesh  point  inside  the  FD  mesh  domain,  i.e. 


E 


^  L~  /  X  /  X  ^  c" 


(42) 


(yc-yr) 


V  ^  (yg-yr)  y  _ 


(yc-ys) 


(yB-yc>  CvB-^r)  (yc-yB>  (yc-y?^ 


^  (yr-y,)  ^ 


(43) 
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C)  Computation  of  Total  Charge 

To  find  the  characteristic  impedance  and  the  phase  velocities  for  a  transmission  line 
having  an  inhomogeneous  medium  requires  calculating  the  capacitances  of  the  structure, 
with  and  without  the  dielectric  regions  [8].  Since  the  capacitance  per  unit  length  is  directly 
related  to  the  charge  per  unit  length  on  the  strips,  the  problem  is  reduced  to  finding  the 
total  charge  per  unit  length  on  the  strips.  The  microstrip  transmission  lines  shown  in  Figs. 
1  to  4  have  two  horizontal  conductors  placed  between  the  two  dielectric  regions.  It  will 
be  assumed  throughout  the  analysis  of  the  problem  that  the  conductors  employed  are 
perfectly  conducting  strips.  It  will  also  be  assumed  that  the  dielectric  materials  are  to  be 
lossless,  isotropic,  and  homogeneous.  Let  Q '  denote  the  total  charge  on  the  f'  conductor. 
If  Gauss’  law  is  applied  to  a  closed  path  /enclosing  the  f  conductor,  the  total  charge  on 
that  conductor  is  then  expressed  as: 


jD’ds  =  ^eE’ds  = 


(44) 


where  D  is  the  displacement  vector  and  ds  is  the  elemental  area.  Substituting  for 
E  =  -VV  in  equation  (44)  yields 


dV  I 

-  VV-n  ds=  -  - dz  dll  = 

s  i  0 


(45) 


The  integrand  in  equation  (45)  is  independent  of  z,  and  hence  the  integral  along  the  z-axis 
is  a  constant.  The  total  charge  per  unit  length  is  then  given  from  equation  (45)  as 
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Q 


me  r 

-  -j.-dt 


g 


dn 


(46) 


The  contour  /  surrounds  the  i‘’'  conductor,  and  it  can  be  explicitly  prescribed  on  the  four 
sides  surrounding  the  conductor;  thus 


Q. 


g 


f  e(x,y) — <fy  +  f  e(x,y) — dc  +  f  e(x,y)—ify  +  f  €(x,y)— <& 
J  dx  J  dy  J  dx  J  dy 

AST  TS  LS 


(47) 


where  the  notation  RS,  TS,  LS,  and  BS  stands  for  the  right,  top,  left,  and  bottom  sides, 
respectively.  The  parameter  g  is  the  length  of  the  transmission  line  along  the  z  direction. 
To  further  illustrate  how  equation  (47)  is  apphed,  consider  the  grid  points  shown  in  Fig. 8. 
For  simplicity,  the  charge  along  the  top  and  bottom  sides  will  only  be  expanded  assuming 
for  a  few  grid  points  on  the  contour  around  the  strip.  Applying  (46)  on  the  right  and  top 
side  yields 


^  (y^-y^)  (ye-y^  j 


(Xj-Xj) 


'(W 

,  (y9-yT> 


(ye-y^  J 


(Xj-Xg) 


((yn-yj  (y>-y,>] 

-  +  - 

^  (yi2-yio>  (yp-y?)  J 


(V  -V 

V  IS  13 
^  Ois-J',,) 


(y,>-y,.) 


■^12  y\()  j 


(^11-^14) 


f2 

(^6-^4) 

- +  - 

(y^-ys) 

2 

t  (»,-».)  (V*.)  J 

2 

(^6-^4)  (^9-^7)  , 

and  on  the  left  and  bottom  sides. 


- +  - 

(yj-^i) 

. 

1 

(Vg-yj) 

^  (Xj-Xj)  (^(5’^4)  J 

*  2 

y  (^6  ; 

(48) 


(49) 


(50) 


Left  Side 


contour ^  ... 


Top  Side 


Fig.  8.  The  FD  mesh  around  the  strip  and  the  labeling  of  the 

points  for  the  application  of  Gauss’s  law  on  the  four  sides 
of  contour  t  in  order  to  compute  the  total  charge  on  the 
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f2 

-  +  - 

(Xc-X.)  + 

((v,-r;i 

-  +  - 

(^g-^s) 

2 

^  Cy3->'i)  (ys-y^  j 

^  (y,-yj)  (y6-y4>^ 

^  f2 

(^11  ■*’ 

- + - - 

2 

^  (yn-yJ  } 

(51) 


where  the  total  charge  per  unit  length  on  the  conductor  is  then  described  as: 


-\qL*  Qts*  QL*  Qbs  . 


»=1,2 


(52) 


D)  Computation  of  Self  and  Mutual  Capacitances,  Impedances,  Phase 
Velocities  and  Coupling  Coefficients 

Following  the  same  procedures  presented  in  [1]  and  [6],  the  characteristics  of  the 
transmission  lines  are  computed.  In  the  previous  section  it  was  shown  how  the  total  charge 
on  each  conductor  may  be  computed.  From  the  total  charge  and  the  voltage  on  each 
conducting  strip,  the  self  and  mutual  capacitances  are  found  for  the  given  transmission 
line.  From  the  total  charges,  the  effective  permitivities  for  all  the  modes  (odd  or  even  for 
the  two  conductor  case)  are  computed,  and  with  them  the  phase  velocities  are  found.  With 
the  phase  velocities,  and  the  self  and  mutual  capacitances  the  characteristic  impedance  of 
the  line  is  computed.  For  the  two  conductor  case,  the  odd  and  even  impedances  are  found, 
and  from  them,  the  characteristic  impedance  for  the  two  conductor  line  is  computed. 

If  a  voltage  V,  on  one  conductor,  excites  a  voltage  Vj  on  the  other  conductor,  then  the 

electric  coupling  which  may  be  written  in  terms  of  the  capacitances  as 
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k 


»12 


=  20  log 


(53) 


But  if  there  is  a  voltage  V2  on  the  second  conductor,  the  one  that  excites  a  voltage  Vj  on 
the  first  conductor  then  the  coupling  equation  is  defined  as 


k 


02] 


=  20  log  JO 


(54) 


For  a  symmetric  microstrip  transmission  line,  Cn  =  C22  ,  C12  -  C21,  and  therefore  kgi2 
=  =  K-  so  the  above  equations  for  the  capacitive  coupling  yield  the  same  result. 


CHAPTER  rV 

NUMERICAL  RESULTS  AND  DISCUSSION 
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A)  Introduction 

This  section  presents  the  results  obtained  from  analyzing  a  series  of  geometries.  Two 
different  FORTRAN  codes  are  developed.  The  first,  generates  the  FD  mesh.  The  second 
solves  for  the  potential  distribution  and  computes  the  line  parameters.  In  this  section,  the 
creation  of  an  input  file  for  the  mesh  generator  is  presented.  Then,  some  results  are  compared 
with  previously  published  data.  Finally,  the  effect  of  the  air  gap  is  shown  for  the  geometries 
shown  in  Figs.  1  to  4. 

B)  Generating  the  FD  Mesh 

A  mesh  generator  program  was  created  to  make  it  easy  to  build  FD  Meshes  for  a 
number  of  geometries  by  writing  an  input  file  that  defines  the  geometry.  To  show  the  way 
such  file  is  created,  the  steps  needed  will  be  followed  for  a  sample  geometry.  The 
geometry  chosen  is  shown  in  Fig.  9.  The  starting  point  for  the  input  data  file  should  be  a 
sketch  of  the  cross  section  of  the  transmission  line  to  be  analyzed  (see  Figs.  1  to  4  and  9.) 
The  dimensions  should  all  be  known  as  well  as  the  permitivities  of  the  materials.  The 
conductors  on  this  geometries  are  assumed  to  be  perfect  electric  conductors  (PECs).  An 
Artificial  Boundary  (ABC)  at  which  the  mesh  will  be  truncated  is  chosen.  The  procedure 
is  the  same  as  that  in  a  SPICE  netlist  where  a  circuit  is  described  by  defining  a  number  of 
nodes  and  then  the  different  elements  between  them  and,  in  this  case,  the  mesh  generator 
input  file  consisting  of  a  series  of  blocks,  their  position,  and  material  is  defined.  There  is 
an  important  rule  when  defining  these  blocks,  that  is,  it  must  have  only  one  kind  of 
dielectric  material,  or  a  PEC.  The  different  blocks  or  cells  on  which  the  cross  section 
of  the  transmission  line  has  been  subdivided  are  shown  in  Fig.  10.  The  Input  Data  File 
follows  the  following  format.  The  first  line  of  the  file  is  a  text  line  describing  the 
geometry,  and  the  next  two  lines  contain  the  number  of  blocks  on  the  x  and  y  directions. 
The  next  set  of  lines  indicate  how  many  points  or  nodes  the  user  wants  on  each  block. 
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Fig.  9.Relation  between  dimensions  of  the  transmision  line 
and  values  of  the  relative  permitivities  used  in  obtaining 
numerical  results,  t/hl  =  0.01,  hl=h2,  s/hl=.5,  W1=W2, 
Wl/hl  =  0.5,  8i  =  2.25,  4.8,  9.6,  82=2.25,4.8,9.6, 
and,  1.0  (no  overlay). 
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Fig.  10.  An  x-y  plane  cut  showing  the  block  structure  to  set  up  the 
input  file  for  the  mesh  generator  program. 
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This  allows  the  program  to  generate,  non  uniform  meshes.  The  input  data  file  appears  as 
follows: 


Line 

#1 

Line 

#2 

Line 

#3 

Line 

#4 

Line 

#5 

Line 

#6 

Line 

#7 

Line 

#8 

Line 

#9 

Line 

#10 

Line 

#11 

Line 

#12 

Line 

#13 

Line 

#14 

Line 

#15 

Line 

#16 

Line 

#17 

Two 

9 

5 

1 

2 

2 

3 

5 

3 

2 

2 

1 

2 

3 

2 

3 

1 


striplines  w/  dielectric  substrate  and  overlay  (Fig,  1.) 

!  Number  of  blocks  on  x  direction 
[Number  of  blocks  on  y  direction 
!No.  of  points  in  the  first  block  in  the  x  direction 


!No.  of  points  in  the  last  block  (#9)  in  the  x  direction 
!No.  of  points  in  the  first  block  in  the  y  direction 


!No.  of  points  in  the  last  block  (#5)  in  the  y  direction 


This  file  must  contain  some  information  about  the  material  and  position  of  each  block 
where  the  next  set  of  lines  are  used  to  represent  material  characteristics.  The  first  line  of 
this  set  is  a  comment  text  line.  This  line  acts  as  a  remainder  of  how  the  next  set  of  lines 
should  be  written.  The  lines  after  the  comment,  contain  a  set  of  integer  and  real  numbers 
which  contain  all  the  information  about  a  given  block.  Each  of  these  lines  defines  a  single 
block.  For  each  line,  the  following  information  must  be  entered:  the  block  index  numbers 
(used  for  identification  of  the  block),  the  lower  left  and  upper  right  comers  of  the  block, 
the  relative  permitivity  of  the  block  (the  Value  0.0  is  used  to  indicate  PEC),  and  the 
potential  for  a  PEC  block  (dielectric  blocks  do  not  need  to  have  a  fixed  potential  specified 
so  0.0  may  be  entered  in  those  cases).  An  example  input  for  Fig.  10.  is  as  follows: 


Line  #18  Block#!  Block#)  xlow  ylow  xup  yup  epsr  volt 
Line  #20  1  1  0.  0.  1.  0.5  1.  0. 

Line  #21  3  2  2.  0.5  3.  2.5  4.4  0. 
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Line  #22  4  3  3.  2.5  4.  3.  0.0  5.  !PEC  strip 


Last  Line  END  !  last  line 

The  blocks  do  not  need  to  be  put  in  any  specific  order  since  the  indexes  of  the  block 
describe  its  location  on  the  cross  section  of  the  transmission  line.  Choosing  some  kind  of 
order,  however,  makes  debugging  easier.  The  output  file  from  this  program  contains  all 
the  information  about  the  geometry  in  a  series  of  points  or  nodes.  To  be  more  precise,  it 
contains:  the  i  andj  indexes  for  each  point;  its  x  and  y  coordinates;  the  type  of  point, 
(general,  boundary,  and  fixed  as  described  in  Chapter  II);  the  permitivities  of  the  four 
regions  around  each  point  (see  Fig.  6.)  or  the  voltage  at  the  point  in  case  of  fixed  potential 
points  (i.e.  points  on  the  strips  or  the  ground  plane). 


C)  Comparison  of  Numerical  Results  with  Previously  Published  Data 

Before  analyzing  any  of  the  geometries  shown  in  Figs.  1  to  4,  it  is  desirable  to  insure 
that  the  developed  code  is  working  properly.  To  accomplish  this,  the  simple  geometry  found 
in  Fig.  3  of  [10]  was  analyzed  for  a  series  of  dimensions.  This  geometry  is  a  two  conductor 
microstrip  line  with  no  overlay,  i.  e.  Cj  is  set  to  1,  free  space,  and  is  set  to  9.6,  AI2  O3.  The 
effective  permitivities  and  the  characteristic  impedance  of  the  line  where  computed  for 
different  dimensions,  special  care  was  taken  in  changing  the  number  of  points  per  block.  By 
doing  so,  the  separation  of  the  nodes  in  the  meshes  was  kept  constant  for  all  the  dimensions. 
The  output  data  obtain  from  the  program  was  compared  with  the  results  from  the  program 
described  in  [10].  Table  I  and  II  contain  these  results. 


31 


Table  I.  Comparison  of  the  results  from  the  FD  code  for  the  characteristic  impedances  for 


Ulc  aiivj 

Dimensions 

Present 

Reference  [10] 

- 

%error 

s/h 

w/h 

Ze 

Zo 

Ze 

z„ 

Ze 

Zo 

0.5 

0.5 

85.97 

46.88 

83.78 

49.85 

2.6 

5.9 

1 

1 

57.14 

42.68 

56.25 

42.85 

1.5 

0.3 

2 

2 

35.64 

32.39 

35.59 

32.08 

0.14 

0.96 

1 

0.5 

78.51 

56.18 

76.92 

57.03 

2 

1.5 

0.5 

1 

61.3 

36.65 

60.22 

37.67 

1.7 

2.7 

1 

0.2 

104.84 

73.61 

103 

77.35 

1.7 

4.8 

Note;  all  the  Impedances  are  expressed  in  Ohms 


Table  II.  Comparison  of  the  results  from  the  FD  code  for  the  effective  permitivities  for  the 
_ even  and  odd  modes  with  the  results  given  by  Smith,  et  al.  in  [10] _ 


Dimensions 

Present 

Reference  [10] 

%error 

s/h 

w/h 

^  eff  even 

^  eff  odd 

^  eff  even 

^  eff  odd 

^  eff  even 

^  eff  odd 

0.5 

0.5 

6.55 

5.15 

6.6 

5.46 

0.75 

5.6 

1 

1 

6.95 

5.62 

6.99 

5.76 

0.57 

2.4 

2 

2 

7.19 

6.39 

7.38 

6.38 

2.6 

0.15 

1 

0.5 

6.6 

5.36 

6.61 

5.56 

0.15 

3.6 

0.5 

1 

6.95 

5.38 

6.99 

5.61 

0.57 

4.1 

1 

0.2 

6.29 

5.18 

6.31 

5.45 

0.32 

4.9 

As  seen  from  these  two  tables  the  results  given  by  the  FD  analysis  are  close  to  those 
reported  in  [10],  The  main  difference  between  the  transmission  line  analyzed  in  [10]  and  the 
one  analyzed  using  the  FD  code  is  the  thickness  of  the  conductors.  While  [10]  assumes 
conductors  with  no  thickness,  the  FD  analysis  assumes  conductors  of  thickness  so  that  t/hl 
=  0.01.  For  all  the  different  dimensions  analyzed,  hi  remained  equal  to  1  unit  and  t  was  set 
to  0.01  units  (see  Fig.  9)  and  for  these  different  cases  only  Wl,  W2  and  S  are  changed.  For 
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the  case  where  s/h  =  W/h  =  0.5,  the  number  of  points  on  the  mesh  was  increased.  Figure  1 1 
shows  the  results  of  these  computations,  in  which  the  convergence  the  effective  permitivity 
and  the  characteristic  impedance  is  very  satisfactory. 

Another  geometry  that  was  analyzed  is  the  one  given  in  Fig.  9  of  [11].  This  geometry 
consists  of  two  microstrips  with  a  groove  etched  between  them.  On  Fig.  11  of  [11],  the 
characteristic  impedance  for  the  even  and  odd  cases  is  plotted  versus  the  depth  of  the  groove 
(hi)  at  94GHz.  From  that  figure  data  is  obtained  for  case  where  hj  =  50.0  pm.  Table  El 
shows  the  results  obtained  using  the  FD  quasi-static  technique  used  in  this  paper  and  the 
results  presented  in  [11]. 


Table  El.  Comparison  of  results  from  the  FD  code  for  the  characteristic  impedances  for 


the  odd  and  even  modes  with  the  results  reported  in  [1 1] 


Present 

Reference  [11] 

Z„ 

43.9Q 

40.5  Q 

Zc 

51.38Q 

54Q 

D)  The  Effect  of  the  Air  Gap  for  the  Large  Overlay  Geometry 
In  the  previous  section  it  was  shown  that  the  output  from  the  program  agrees  with 
previously  published  data.  Now  the  effect  of  the  air  gap  was  studied.  The  geometry  with  the 
large  dielectric  overlay  (see  Figs.  1  and  3)  was  analyzed  for  the  dimensions  given  in  Fig.  9. 
Different  values  of  permitivities  for  both  the  substrate  and  the  overlay  were  used.  The  code 
must  first  solve  for  the  potential  distribution  on  the  mesh  before  computing  all  the  line 
parameters.  The  different  potential  distributions  for  the  odd  mode  for  the  geometries  shown 
in  Figs.  1  and  3  are  shown  on  Fig.  12a.  The  dimensions  are  those  given  in  Fig.  9,  the 
permitivities  are  set  so  that  —  9.6  and  =  4.8.  It  is  clear  from  Fig.  12a  that  the  potential 
on  the  strips  is  fixed  to  1  volt  or  -1  volt.  The  ground  plane  is  clearly  seen  set  to  zero  potential. 
It  is  also  clear  that  the  chosen  ABC  is  working  properly.  There  is  not  a  clear  difference 
between  the  potential  distribution  for  the  geometry  with  an  air  gap  present  (Fig.  1)  and  the 
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Fig.  11.  Convergence  plot  with  no  overlay,  8,1  =  9.6, 
8,2  =  1,  W1=W2,  and  S/hl  =  Wl/hl  =  0.5. 
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Fig.  12a.  Potential  distribution  on  the  FD  mesh  for 

=  9.6,  8,2  =  4.8,  ^1=  -1  Volt,  and  ¥^=1  Volt. 
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one  without  one  (Fig.  3).  The  potential  distributions  for  the  even  mode  are  shown  on  Fig  12b. 
The  strips  with  their  potential  fixed  at  1  volt  are  easily  identified.  The  ground  plane  set  at 
zero  volts  is  also  clear.  Right  above  the  strips  the  interface  between  the  overlay  and  the  free 
space  on  top  of  it  can  be  seen.  The  difference  between  the  case  with  an  air  gap  and  without 
air  gap  is  a  little  more  clear  than  in  the  odd  mode.  The  potential  at  a  point  between  the  strips 
is  lower  when  the  air  gap  is  present.  For  this  same  case  the  number  of  nodes  on  the  mesh  was 
increased  to  check  the  stability  of  the  code.  Fig.  13  shows  the  convergence  of  the  capacitive 
coupling  and  the  self  and  mutual  capacitances,  that  are  equal  for  a  symmetric  geometry.  Now 
for  the  geometries  shown  on  Fig.  1  and  3,  the  permitivity  of  the  substrate  was  left  constant 
and  the  permitivity  of  the  overlay  was  changed.  Figs.  14a  to  14c  show  the  capacitive 
coupling  and  the  mutual  and  self  capacitances  as  a  function  of  the  overlay  permitivity.  For 
all  the  cases,  the  presence  of  the  air  gap  reduces  the  capacitive  coupling  since  the  capacitance 
between  the  strips  is  reduced.  The  electric  field  distribution  in  the  mesh  is  also  obtained. 
Figures  15a  to  15d  show  the  vector  field  distribution  for  the  even  mode,  with  and  without  the 
air  gap.  Looking  at  the  entire  mesh,  it  does  not  appear  to  be  any  difference  on  the  distribution 
of  the  electric  field.  Figures  15b  and  15d  show  the  area  between  the  conductors,  here  the 
effect  of  the  air  gap  is  clearly  seen.  Figures  16a  to  16d  show  the  distribution  for  the  odd 
mode.  Again  no  effect  is  seen  when  the  distribution  on  the  whole  mesh  is  studied.  The  details 
of  the  area  between  the  conductors  (Figs.  16b  and  16c)  show  a  small  effect  on  the  electric 
field,  but  not  as  big  as  the  one  for  the  even  mode  (Figs.  15b  and  15d). 


E)  The  Small  Overlay  Geometry 

The  small  overlay  geometry  shown  on  Figs.  2  and  4  has  also  been  analyzed.  The 
dimensions  used  were  those  of  a  real  microstrip  transmission  line  that  was  used  as  a  probe 
to  measure  the  relative  permitivity  of  unknown  materials.  The  probe  has  the  following 
dimensions;  the  width  of  both  conductors  is  4mm,  separated  by  2.5mm,  while  their  thickness 
was  calculated  from  the  manufacturers  specification  of  the  amount  of  copper  per  square  inch. 
The  calculations  yield  0.071 12mm,  and  the  same  thickness  was  use  to  represent  the  ground 
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Fig.  13.  Convergence  of  k  and  the  self  and  mutual  capacitance 
for  the  large  overlay  geometry  =  9.6,  8^2  =  4.8 
and  with  an  air  gap. 


Fig.  14a.  Coupling  coefficient,  self  and  mutual  capacitances  for 
different  overlays,  with  and  without  air  gap  with  e„  -  y.O. 
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Fig.  14c.  Coupling  coefficient,  self  and  mutual  capacitances  for 
different  overlays,  with  and  without  air  gap  with  =  2.25. 
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Fig  15  a.  Vector  map  of  the  electric  field,  large  overlay 
geometry,  even  mode,  air  gap  present. 


Fig  15  b.  Detailed  vector  map  of  the  electric  field  between  the 
conductors,  large  overlay  geometry,  even  mode,  air 
gap  present. 
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Fig  15  c.  Vector  map  of  the  electric  field,  large  overlay 
geometry,  even  mode,  no  air  gap  present. 


Fig  15  d.  Detailed  vector  map  of  the  electric  field  between  the 
conductors,  large  overlay,  even  mode,  no  air 
gap  present. 
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Fig  16  a.  Vector  map  of  the  electric  field,  large  overlay 
geometry,  odd  mode,  air  gap  present. 


Fig  16  b.  Detailed  vector  map  of  the  electric  field  between  the 
conductors,  large  overlay  geometry,  odd  mode,  air 
gap  present. 
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Fig  16  c.  Vector  map  of  the  electric  field,  large  overlay 
geometry,  odd  mode,  no  air  gap  present. 


Fig  16  d.  Detailed  vector  map  of  the  electric  field  between  the 
conductors,  large  overlay  geometry,  odd  mode,  no  air 
gap  present. 
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plane.  The  height  of  the  substrate  was  found  to  be  1.651  mm,  and  the  overlay  was  set  to 
0.8mm.  The  substrate  extends  5mm  from  the  outer  edges  of  the  conductors,  while  the  overlay 
extends  only  2.5mm.  The  relative  permitivity  of  the  substrate  was  kept  at  4.8  while  the 
relative  permitivity  of  the  overlay  was  changed.  Like  in  the  previous  cases  for  the  large 
overlay,  the  presence  of  the  air  gap  did  reduce  the  coupling.  Fig.  17  show  the  results  obtained 
for  the  self  and  mutual  capacitances  and  the  coupling  coefficient  with  and  without  the  air 
gap..  Table  HI  shows  the  results  for  the  case  when  the  relative  permitivity  of  the  overlay  is 
9.6. 


Table  IV.  Results  for  the  small  overlay  geometry. 


Cn 

Ci2 

K 

No  air  gap 

282.67 

35.93 

-18.97dB 

Air  Gap 

297.1 

30.3 

-20.39dB 

Figure  18a  shows  the  electric  field  around  the  conductors  for  the  even  mode.  The  effect  of 
the  air  gap  is  clearly  seen  and  it  appears  to  be  stronger  for  the  even  mode  than  for  the  odd 
mode  shown  in  Fig.  18b. 

F)  Coupling  on  an  Asymnnietric  Transmission  Line. 

Figure  19  shows  two  different  geometries  that  were  also  analyzed.  The  top  one  has 
the  conducting  strips  embedded  in  the  dielectric  substrate.  The  bottom  one  has  the  conductors 
resting  above  the  substrate.  Furthermore,  the  substrate  contains  different  types  of  dielectric 
materials.  The  dimensions  of  the  geometry  were  the  following:  The  height  of  the  substrate  (h) 
was  0.3mm,  the  width  (w)  of  the  conductors  was  0. 1mm,  the  separation  between  the  lines  (s) 
was  0.075mm,  and  the  thickness  of  the  conductors  (t)  was  set  to  0.001mm.  The  substrate  was 
composed  of  the  following  materials:  polystyrene,  silicon,  and  alumina  with  dielectric 
constants  of  2.54,  1 1.9  and  9.6,  respectively.  The  values  where  taken  from  Appendix  G  of 
[12].  Table  V  shows  a  sample  of  numerical  results  for  both  geometries  shown  in  Fig.  19. 
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Table  V.  Results  from  the  FD  code  for  the  geometries  shown  in  Fig.  19. 


Case 

Numerical  results  based  on  the  FD  code 

Erl 

er2 

^r3 

Embedded  lines 

Kn 

Kii 

^eff  even 

^efifodd 

2.54 

9.6 

9.6 

no 

-18.96dB 

-18.96dB 

5.14 

3.77 

2.54 

11.9 

11.9 

no 

-20.15dB 

-20.15dB 

6.11 

4.36 

2.54 

9.6 

11.9 

no 

-18.87dB 

-20.22dB 

5.14 

3.78 

2.54 

9.6 

9.6 

yes 

-18.33dB 

-18.33dB 

5.16 

3.86 

2.54 

11.9 

11.9 

yes 

-19.52dB 

-19.52dB 

6.14 

4.44 

2.54 

9.6 

11.9 

yes 

-18.25dB 

-19.59dB 

5.16 

3.86 

These  results  show  that  the  presence  of  air  between  the  conductors  reduces  the  coupling 
between  the  lines  by  approximately  0.5dB.  The  reduction  is  not  significant  compared  to  the 
previous  cases  where  an  overlay  was  present.  The  increase  in  the  dielectric  constant  of  the 
substrate  under  the  conductor  increases  the  self  capacitance  on  the  line.  This  increase  in  the 
self  capacitance  translates  in  a  reduction  of  the  coupling  between  the  lines. 
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£2 

Fig.  17.  Coupling  coefficient,  self  and  mutual  capacitances  for 
different  small  overlays,  with  and  without  air  gap. 
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Fig.  18  a.  Detailed  vector  map  of  the  electric  field  around  the 
conductors  for  the  small  overlay  geometry  with 
and  without  the  air  gap,  even  mode. 
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Fig.  18  b.  Detailed  vector  map  of  the  electric  field  around  the 
conductors  for  the  small  overlay  geometry  with 
and  without  the  air  gap,  odd  mode. 
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Fig.  19.  Two  geometries  with  substrates  made  out  of  different 

dielectrics.  One  has  conductors  embedded  in  the  substrate 
and  the  other  has  conductors  lying  on  the  substrate. 
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CHAPTER  V 
CONCLUSION 

The  FORTRAN  code  developed  can  analyze  a  large  variety  of  transmission  line  geometries. 
The  results  that  have  been  compared  with  previously  published  were  very  satisfactory.  The  program 
has  shown  that  it  can  deal  with  complicated  geometries  containing  different  dielectric  materials  and 
that  these  geometries  may  be  asjmimetric.  The  effects  of  the  presence  of  the  air  gap  on  the  electric 
field  distribution  have  been  shown.  These  changes  lead  to  a  reduction  of  the  capacitive  coupling 
between  the  lines.  However,  the  coupling  may  be  increased  when  the  relative  permitivity  of  the 
overlay  is  large.  When  building  couplers  it  may  be  desirable  to  eliminate  the  air  gap  by  embedding 
the  lines  in  the  substrate  or  the  overlay.  It  is  also  shown  that  susbtrates  with  high  dielectric  constants 
help  reduce  the  coupling. 
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